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We present an STZ-based analysis of numerical simulations by Haxton and Liu (HL). The ex- 
tensive HL data sharply test the basic assumptions of the STZ theory, especially the central role 
played by the effective disorder temperature as a dynamical state variable. We find that the theory 
survives these tests, and that the HL data provide important and interesting constraints on some 
of its specific ingredients. Our most surprising conclusion is that, when driven at various constant 
shear rates in the low-temperature glassy state, the HL system exhibits a classic glass transition, 
including super- Arrhenius behavior, as a function of the effective temperature. 



I. INTRODUCTION: BASIC HYPOTHESES 

In a remarkable recent report, T. Haxton and A. Liu 
(HL) [l[ have described molecular-dynamics simulations 
of a simple, sheared, two-dimensional, glass-forming ma- 
terial over three decades of steady-state strain rates 7, 
and for bath temperatures T ranging from about one 
tenth of the glass transition temperature To to about 
twice To. Most importantly, by measuring pressure fluc- 
tuations, HL have determined values of the effective dis- 
order temperature T e g- for each value of 7, T, and the 
shear stress s. This extensive data set tests the applica- 
bility of any theory of amorphous plasticity such as the 
shear-transformation-zone (STZ) theory discussed here, 
and also probes the limits of validity of the effective- 
temperature concept. 0, H, 0, H, IE 0| We find that the 
main features of the STZ theory presented in earlier pa- 
pers [1, H, [l(| [LH nicely survive these tests. Moreover, 
several of the theory's specific ingredients, which had not 
been tightly constrained by earlier experiments or simu- 
lations, can be refined and extended in light of the HL 
data. Our most surprising conclusion is that, for T < Tq, 
the relation between 7 and T e g- is a direct analog of the 
relation between the a relaxation rate and the bath tem- 
perature T near a conventional glass transition. 

It is well-known (e.g. see [12|]) that, for unstressed 
glass-forming materials in a range of temperatures above 
some Ta, the viscosity r\ is - at least to a first approxi- 
mation - a simple activated process with a temperature- 
independent energy barrier. Below Ta, T) depends super- 
Arrheniusly on T, diverging as T decreases toward the 
glass transition temperature (or Kauzmann temperature) 
Tq. Figure [lja) shows that the HL data for the viscos- 
ity in the limit of vanishing strain rate can be fit by a 
function of the form 



our data analysis in Section IHI 



i] = 770(0) exp 
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a{6) 



(1.1) 



where 9 = T/Tz and kgTz is an Arrhenius activation 
energy that we interpret as the energy of formation for an 
STZ. a{9) is a model-dependent function, to be specified 
later, that diverges as 9 — > 9q = Tq/Tz and vanishes for 
9 > 9a = Ta/Tz- The prefactor 770 (9) is a relatively 
slowly varying function of 9 that appears prominently in 
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FIG. 1: (a) Logarithm of the HL viscosity (open triangles) as 
a function of the inverse temperature 1/6 — Tz/T. The solid 
line is the STZ-theory fit to these data points. The dashed ex- 
tension of this line indicates a region where the temperatures 
are so high that the theory becomes incorrect, (b) Logarithm 
of the HL strain rate as a function of the inverse effective 
temperature 1/x at temperatures below the glass transition 
temperature 80 = 0.20. The data points correspond to tem- 
peratures — 0.022 (stars), 0.044 (squares) and 0.11 (trian- 
gles). The solid line is the theoretical fit. See Section ITTT1 for 
details of the theoretical analysis and an evaluation of Tz ■ 

Our key hypothesis is that, below the glass transi- 
tion temperature 9q, there is a relationship analogous to 
Eq. ljl.ip between the steady-state effective temperature 
and the strain rate 7. In analogy to the definition of 9, 
and consistent with earlier notation [9|, llfj . Illj , we write 
X = T e g/Tz- We also define q = tq 7, where To is the fun- 
damental time scale that appears in the dynamical STZ 
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equations of motion. Thus, q is a dimensionless rate of 
molecular rearrangements or, equivalently, a stirring rate. 
The HL results are expressed in terms of a molecular time 
scale, to i which is derived from their molecular interac- 
tions and is not necessarily the same as To - although the 
two times should have about the same order of magni- 
tude. We find it convenient to set tq = to = 1, and to 
incorporate any difference between these time scales into 
the parameter eo defined below in Eq. (|2.5p . For clarity, 
however, we occasionally retain explicit factors r . 

We propose that the rate q be interpreted as a purely 
kinetic quantity, independent of the specific mechanism 
by which the rearrangements occur. We then define x(q) 
to be the dimensionless, steady-state, effective temper- 
ature at stirring rate q. To a first approximation, we 
assume that xil) ls independent of 9 for 9 < 9q. As 
we shall see later, the temperature dependence of x(l) 
becomes more interesting for larger 9. 

In the limit g C 1, the number of rearrangements de- 
pends only on the extent of the deformation and not 
its duration; thus the condition for steady-state equilib- 
rium is that, with high probability, each molecule has 
changed its neighbors at least a few times. It follows that 
x(q — > 0) = xo is a nonzero constant, also independent of 
the rearrangement mechanism. We then postulate that 
the inverse function q(x) is the rate at which rearrange- 
ments occur when the system is driven so that its steady- 
state effective temperature is x- Then, in analogy with 
Eq. (|l.l|) , we write 
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(1.2) 



where qo and A are constants to be determined, and 
a eff(x) is a super- Arrhenius function that diverges at the 
effective Kauzmann temperature xo and vanishes above 
some xa- 

This hypothesis is supported by the HL data. Fig.QJb) 
is an Arrhenius plot of the HL results for the strain rate 
versus the effective temperature for three values of 9, all 
almost certainly below 9$, and the solid curve is a fit to 
those points using Eq. (jl.2|) . There are large uncertain- 
ties in the data at small strain rates (large 1/x), but the 
existence of Arrhenius and super- Arrhenius regions is ap- 
parent. Note that the Arrhenius region, i.e. the section 
of the curve that is linear in 1/x, seems to extend all the 
way to infinite x- Ordinarily we do not think of Arrhe- 
nius formulas remaining valid out to high temperatures. 
The difference here may be that the effective tempera- 
ture pertains only to a fraction of the degrees of freedom 
in this system; so, perhaps in analogy to magnetic spin 
temperatures, the maximally disordered state at infinite 
X is realistic. 

One of the most interesting features of Eq. (|1.2[) is that 
X diverges at a finite strain rate, qo = 0.08 as seen from 
the extrapolation shown in Fig. QJb) . For molecular sys- 
tems where r is of order femtoseconds, 7 = 0.08/r is 
an impossibly large rate; but it is accessible in numerical 
simulations and perhaps also in granular systems or the 



like. The system must liquify in some sense when driven 
faster than qo; and the STZ theory must break down near 
that point. 



II. STZ THEORY 
A. General STZ Equations 

In order to use Eq. (|1.2p to compute stresses, strain 
rates and effective temperatures for 9 both below and 
above 9q, we need a complete set of dynamical equa- 
tions of the kind provided by the STZ theory. In this 
Section, we summarize the version of that theory that 
was presented in [ll[ , specializing to steady-state motion 
but otherwise retaining the generality needed to describe 
plasticity over the wide range of driving conditions stud- 
ied by HL. We need one structural modification of this 
theory, which we explain following Eq. (|2.8p below; and 
our analysis requires model dependent rate factors that 
differ from those used in earlier papers. We start by writ- 
ing the STZ equations in a form determined by the basic 
structure of the theory. Then, in the following subsection, 
we specify our choice of the model dependent ingredients 
of these equations. 

In keeping with the emphasis of HL and the results of 
recent STZ analyses such as d, [loj], we write this the- 
ory in a form in which the effective temperature x is the 
single, dynamically relevant, internal state variable. Ac- 
cordingly, we need just one equation for the dimensionless 
strain rate q as a function of the shear stress s and the ef- 
fective temperature x? supplemented by the steady-state 
version of an equation of motion for x- 

The equation for the strain rate has the form 



7 T = 2 e-Vxf&e), 



(2.1) 



where f(s,9) is a model-dependent function (to be spec- 
ified shortly), and exp(— 1/x) is proportional to the 
STZ density. This clean separation between the stress- 
dependent and x-dependent parts of the plastic strain 
rate is an important characteristic feature of the STZ 
theory that has been emphasized recently by Shi et. al. 

The equation for x is 



- 1/x T(s,9) 



1 - 



X 



= K 



(x)p(0) (J 
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(2.2) 



This is essentially the steady-state version of Eq.(3.5) 
in [HI; but here we have written it in a form that is 
more appropriate for situations in which x(q) can become 
large. 

According to Eq. (|2.2[) . x finds its steady-state value 
when the rate at which it is being driven mechanically to- 
ward its kinetically optimal value x(q), on the left-hand 
side, is balanced by the rate at which x relaxes toward the 
bath temperature on the right. T(s,9), another model- 
dependent function, is proportional to the rate at which 
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the work done by the driving force is dissipated as con- 
figurational entropy; thus r(s, 9) vanishes at s — 0. The 
function 



p{0) 



^ J exp [-a(0)] for 9 > 9 







for 9 < 9o 



(2.3) 



is the super- Arrhenius part of the a relaxation rate. Be- 
cause p{9) vanishes for 9 < do, Eq. (|2.2[) implies that 
the low-temperature, steady-state value of x is x(?) as 
anticipated in Eq. (|1.2p . The function K(x) appearing 
in Eq. (|2.2|) is proportional to the density of sites at 
which thermal equilibration events - as opposed to shear- 
transformations - take place. We assume here, as in pH . 
that 



k e 



-P/x 



(2.4) 



with (3 = 1, which means that both kinds of sites are 
comparably populated. The parameter k remains to be 
determined from the data. 

In the general STZ equations as developed in [l(J EH , 
the function f(s,9) in Eq. (|2.ip has the form 



f(s,9) = e C(s) \T(s)-m(s) 
where ep is a dimensionless constant, 



C(s) 



R(s) + R(-s) 



T{s) 



R(s)-R(-s) 
R{s) + R{-s) : 



(2.6) 



(2.7) 



and R(s)/tq is the rate at which individual STZ's un- 
dergo forward shear transitions driven by the stress s. 

The function m(s) is the orientational variable that 
emerges in the role of a back stress in the STZ theory. It 
satisfies the steady-state relation: 

2C(s) T(s)-m(s) l-m(s)s/s =m(s)p(9). (2.8) 

Here, and in the equation for T(s, 9) that follows, we de- 
part slightly from the equations derived in In that 
paper, the yield stress s y was assumed to be approxi- 
mately equal to a temperature-independent constant p, 
and the stress was everywhere written in units of that 
constant. Here, on the other hand, it is obvious - most 
visibly in the HL graphs of stress versus strain rate at low- 
temperatures - that s y must be a function of 9. Exami- 
nation of the derivation of the STZ equations in [TTf| in- 
dicates that the factor of proportionality between T(s, 9) 
and the energy dissipation rate per STZ, which was set 
equal to p, there primarily for dimensional reasons, should 
have been a possibly temperature-dependent stress that 
we now call sq. The result of this change is that a fac- 
tor s, which multiplies the strain rate in the expression 
for the rate at which plastic work is done on the system, 
is replaced by s/sq, and the explicit s that appears in 



Eq. (|2.8[) now becomes s/sq. This change in the analysis 
removes an ambiguity in the older STZ theories; e.g. see 
[3] Details will be published elsewhere. [HI 

The significance of the quantity sq is that below the 
glass transition temperature it is closely related to the 
dynamical yield stress s y . For 9 < 9q and p(9) — 0, s y is 
the value of s for which the two factors on the left-hand 
side of Eq. (|2 . 8[) vanish simultaneously; therefore it is the 
solution of the equation 



s y T(s y 



s - 



(2.9) 



Because T(s) < 1 in all cases and is approximately equal 
to unity at low temperatures, s y = sq. If sq is a function 
of 9 in Eq. (|2.9p . then s y also is temperature dependent. 
Although there is no yield stress for 9 > 9o, f(s,9) still 
changes abruptly at s = so for 9 not too much larger than 
#o- This implies that at stresses close to ,s , the stress vs. 
strain rate curve exhibits an abrupt change in slope. 

The STZ formula for the dissipation rate F(s, 9) ap- 
pearing in Eq. (|2.2[) is 



(2.5) r( s ,0) 



2C(s) 


T(s) — m 




s/s Q - £(m) 


+ m£(m) p{9) 




. — m£(m) 



where m = m(s) as given by Eq. (|2.8[) . and 
so£( m ) = T _1 (m) 



(2.10) 



(2.11) 



is the inverse of the function T(s); that is, T[sg t;(m)] = 
to. It is not difficult to show that eoT(s,9) = 2sf(s,9) 
for 9 < 9q , confirming that the energy dissipated is equal 
to the work done on the system when configurational 
relaxation cannot occur via thermal fluctuations alone. 
Above the glass transition, on the other hand, the exter- 
nal work is not the same as the dissipation, and the use 
of Eq. (|2.10p is necessary for accurate calculations. 



B. Model-Dependent Ingredients 

We turn now to the STZ forward transition rate R(s). 
Here we depart from the choices made in earlier papers 
[ToL [lil ] , where we were concerned with experiments and 
simulations performed only at relatively small driving 
forces, and which therefore did not sharply constrain our 
choice of R(s). HL have provided a data set that goes up 
to stresses more than ten times the yield stress, at tem- 
peratures both well above and well below the glass tran- 
sition, and thus have made it necessary for us to consider 
this rate factor in a broader context than before. 

We need a function R(s) that interpolates between 
Eyring-like behavior at very small stresses and power- 
law growth at very large ones. In choosing such an i?(s), 
we find that we must depart from our earlier, purely 
phenomenological procedure of choosing simple functions 
with very few adjustable parameters, in the expectation 
that whatever data we had available would not justify 
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additional theoretical complications. Here we are facing 
a very different situation. In order to interpret the HL 
data, we find it better to start with a physically moti- 
vated model containing an uncomfortably large number 
of parameters. Determining the values of these param- 
eters is made easier in places by the fact that we can 
fit some parts of the data independently of other parts, 
and then see how these initially separate pieces of the 
puzzle fit together, for example, in connecting our pre- 
dictions of behaviors at very small strain rates to those at 
very large ones. Another advantage of this more nearly 
first-principles procedure is that we can guess the magni- 
tudes of many parameters on physical grounds, and also 
know the limits of validity of physics-based approxima- 
tions. For example, we know that the Eyring-like formula 
that we use for the STZ transition rate breaks down at 
the highest temperatures simulated by HL; and thus we 
know to assign less weight to our theory at those tem- 
peratures in fitting the HL viscosity curve. 
Our proposed R(s) is 



R(s) 



exp 



1 + I — 

si 



i/2 



(2.12) 



The first factor on the right-hand side of Eq. (|2.12p is the 
Eyring rate in a form similar to that used in @] , where the 
exponential function of s/p, causes the rate to saturate 
at large s. Here, 9e is the height of the Eyring activation 
barrier in units of Tz- We expect Oe < 1, because the 
barrier opposing STZ transitions from one state to the 
other should be less than the barrier opposing creation of 
new STZ's. The parameter jl is related to the curvature 
of this barrier at its peak. It appears here ostensibly in 
the same place where Ji appeared in pd| . but now it has 
no direct connection to the yield stress. In the limit of 
small stress, we recover an Eyring rate factor in the form 



R{s) ps exp 



(1 - s/p) 



therefore 



and 



C (s) «e^cosh(^) 



(2.13) 



(2.14) 



stress sq(9) so that sq(0) — 1 and measure stresses s in 
units of this zero-temperature yield stress. The crossover 
stress si determines when the system departs from Eye- 
ring stress-activated behavior. 

The exponent n in Eq. (|2.12p is especially interest- 
ing. In strongly dissipative systems, we expect n = 1, 
indicating that the STZ transition rate is controlled 
at large stresses by some linear dissipative mechanism 
such as friction or viscosity. In the opposite limit, 
where rates are controlled by collisions between primar- 
ily hard-core objects and the detailed molecular interac- 
tions are relatively unimportant, we expect to find Bag- 
nold scaling. [l6| That is, if there is no natural energy or 
stress scale in the problem, then dimensional analysis re- 
quires that the stress be proportional to the square of a 
rate. Since R(s)/tq is the only available quantity with 
dimensions of inverse time, we must have R(s) ~ s 1 / 2 
and n = 1/2. It appears that the latter case is more 
consistent with the HL viscosity data. 

Other model-dependent ingredients that we must spec- 
ify are the super-Arrhenius function a(9) and its ana- 
log a e g(x)- Here we adopt formulas motivated by the 
excitation-chain theory of the glass transition [lj, [l8| 
proposed recently by one of us. For present purposes, 
these formulas need not be interpreted as anything more 
than phenomenological fits to experimental data. We 
write 







exp 



-b- 



(2.17) 



where 6 = 3 produces a reasonably accurate description 
of the transition from an inverse power-law divergence 
near 9 Q to simple Arrhenius behavior above 9a in analy- 
ses of experimental data. [l5j The excitation-chain theory 
predicts p — 1 (Vogel-Fulcher) in three dimensions and 
p = 2 in two, but the HL data is not accurate enough to 
distinguish between these possibilities. Since HL fit their 
two dimensional viscosity data with a Vogel-Fulcher func- 
tion, we have chosen to be conservative and use p = 1 in 
our analyses. However, we point out in Section HTT1 that 
p = 2 remains an interesting possibility. 

The effective-temperature analog of Eq. (|2.17p is 



Xi 



X-Xo 



exp 



X - Xo 



XA - Xo 



(2.18) 



T(s) ps tanh 



IE S 

9~P~ 



(2.15) 



The second factor on the right-hand side of Eq. (|2.12|) 
converts the saturated Eyring function at large s to a 
power law: 



C(*)«(-) ; 



(2.16) 



Here, S\ is a temperature independent stress scale. With- 
out loss of generality, we normalize the nominal yield 



III. COMPARISONS BETWEEN STZ THEORY 
AND THE HL SIMULATIONS 

A. Evaluating the Parameters 

We turn now to a detailed numerical comparison be- 
tween these theoretical predictions and the HL data. We 
start by deducing values of the parameters Tz and A in 
Eq. (fL~2")) and the exponent n in Eqs. (j2~T2"]) and (|2~TB)) . us- 
ing primarily the data in the extreme Arrhenius region 
at low temperatures 9 and large values of x- 
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FIG. 2: Effective temperature \ as a function of log 10 (g) for 
HL results both above (open symbols) and below (nlled-in 
symbols) the glass transition. The smooth curves are the 
predictions of the STZ theory at the same temperatures 9. 

Figure QJb) shows the HL data for — log 10 (q) as a 
function of 1/x, for the three lowest HL temperatures, 
9 = 0.022, 0.044, and 0.11, all of which are comfort- 
ably below the glass transition temperature 9o — 0.20. 
Therefore, p(6) = and x = X- From Eq. p.2p and 
the observed slope of this graph in the Arrhenius region, 
i.e. from the data shown in HL Fig. lb, we find that 
AT Z = .0068 (in HL units). 

A second piece of information is obtained by looking 
at the HL results for stress as a function of 1/x, shown 
here in Fig.[5j Again, we look only in the large-x, small 9 
limit. Equations (|1.2[) and (|2.1[) . plus our large- s estimate 
for the rate factor, combine to give 

q^q e- A ^^2e e- 1 ^(-£\ (3.1) 

from which we find 

« - In (£L) _ ±=±. (3.2) 

n \ 2 e J nx 

The observed slope in Fig. [5J i.e. in the original HL figure 
in which In (s) is plotted as a function of 1/T e g-, tells us 
that (A — X)Tz/n= .004. 

A third relation, equivalent to the first two but useful 
for data analysis, is 

b w s ''4 h (^) + s(^ i ) h(,) ' ( "» 

The observed asymptotic slope of the graph of ln(s) ver- 
sus ln(q) shown in Fig. [6] is (A — l)/nA = 2/3. If we 
assumed n = 1, we would find A = 3 and Tz = 0.00225, 
a value of Tz that is too small to be consistent with the 
HL viscosity data. However, n — 1/2 (the Bagnold pre- 
diction) implies that A = 1.5 and Tz — 0.0045, which 
fits the viscosity quite well. 

We then use the last two graphs to evaluate the q- 
independent terms on the right-hand sides of Eqs. (|3.2|) 



FIG. 3: The same data as shown in Fig[2]but in the form of an 
Arrhenius plot of log 10 (g) as a function of 1/x, for HL results 
both above (open symbols) and below (filled-in symbols) the 
glass transition. The smooth curves are the predictions of the 
STZ theory at the same temperatures 8. 

and (|3.3[) . with q = .08. In this way we check for con- 
sistency of the estimates we have obtained from a noisy 
data set, and find a relation between s± and eo: 

si 2.3 x 10 3 4- ( 3 - 4 ) 

Several other parameters can be determined by direct 
examination of the HL data. HL measure their stresses a 
in units related to their interatomic forces. Because our 
values of s are expressed in units of the zero-temperature 
yield stress, the ratio of a to s is approximately the value 
of a at the lowest strain rate and the lowest temperature 
shown in HL Fig. la, i.e. approximately 0.001 in HL stress 
units. Thus, the stress-conversion relation is 

a = a HL s- &hl — 0.001. (3.5) 

Similarly, we deduce from the small-g limits of the stress 
in HL Fig. la (shown here in Fig. that, for 9 < 6q, 
Sy — so is roughly a linearly decreasing function of 9. 
There is no comparably systematic way to evaluate sq(9) 
for 9 > 9q; but we see no structure in these curves that 
might indicate a further decrease of the stress scale Sq(9) 
at temperatures above #o- Accordingly, we choose 

So{0) = \l-c fortf>fl , ^ 

where c = 0.6. 

The parameter eo, in our notation, is proportional to 
the number density of STZ's in the limit of infinite x> 
modified by a time-scale conversion factor of order unity. 
We have found no independent estimate of eo compara- 
ble to the preceding estimates of Tz, A, n, etc.; but we 
expect it to be small if the STZ theory is to remain valid 
at large values of x- Eq. (13. 4[) provides an additional 
constraint. As discussed in the next section, the stress 
crossover behavior for the full set of HL data provides an 
estimate of Si = 0.08, and this implies that eo = .006. 
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Our next step is to use the parameters determined so 
far to compare the HL viscosity shown in Fig. QJ, with 
the STZ predictions and thereby evaluate the parameters 
that occur in Eqs. flOp , (j2~T51) . and (j?T7j) . Figure H 
shows the HL data for log(s/q) as a function of log(g) at 
temperatures above and below the glass transition. By 
definition, the HL viscosity is 7/ = uhl s /j in the limit 
s — > 0, and the HL viscosity data in Fig. QJi are based 
on extrapolation of the data in Fig. [4] to s = using the 
Cross equation. 
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FIG. 4: Logarithm of the stress s divided by q as a function 
of log (9) for HL results at temperatures both above (open 
symbols) and below (filled-in symbols) the glass transition. 
Lhe curves are the predictions of the SLZ theory at the same 
temperatures. For temperatures below the glass transition, 
the limit as s — > of s/q is the Newtonian viscosity, rj. 



To evaluate 77, we use Eqs.lpTIj). (I23j) . and (pT8j) with 
X = 9, and, after some algebra, find that 



<?HL 

eoT'(0) 

<JHL /i 

£0 Se 



1 



m 

2C(0) 

\w* at9 



(3.7) 



The second version of this equation makes explicit use 
of Eqs. (|2. 14[) and (|2.15|) . Obviously there are too many 
parameters here for us to obtain more than a very rough 
estimate of their range of values. None of the six HL 
data points for rj are close enough to the transition tem- 
perature 9o to provide an accurate picture of the super- 
Arrhenius singularity. To make matters worse, the high- 
est reported temperature is approximately equal to Tz, 
well above any reasonable value for 9e and thus far too 
high for the Eyring formula to be valid. Nevertheless, we 
find some interesting information by exploring the ways 
in which we can fit the theory to the data. 

Our procedure has been simply to use Eq. (|3.7[) in con- 
junction with Eq. ljl.ip to plot the logarithm of rj(6) and 
compare the results with the HL data points. To start, 
we accept the values Tz = 0.0045 (in HL units), obtained 
from the large-stress behavior, and #0 — 0.20 estimated 
by HL. We also fix p = 1 and b = 3 in Eq. (|2~T7|) . Then 
our best-fit parameters are 9a — 0.69, 0\ = 0.69, p, = 0.3, 
and Be — 0.4. Our fit to the HL viscosity measurements 



is shown in Fig. [lja). The dashed part of that curve 
indicates the region where, given the above value of 9e, 
we know that the Eyring formula in Eq. (|2.12p cannot be 
valid. 

In contrast to the uncertainties that complicate the 
theoretical fit to the HL viscosity data, the analogous 
relation between the strain rate and the effective tem- 
perature shown in Fig. [Tf b) seems relatively easy to un- 
derstand. Whereas the viscosity is an intrinsically dy- 
namic quantity, involving material-specific relations be- 
tween stress and plastic response, the function q(x) ap- 
pears to be, as we have postulated, a purely kinetic 
relationship. The solid curve shown in Fig. [TJb) is a 
fit to Eq. (|1.2[) with the following parameters defined in 
Eq. (EUD; xo = 0.20, X a = 0.30, X i = 0.26, p = 1, and 
b = 3. We emphasize that there are still substantial un- 
certainties in these numbers. Within these uncertainities, 
it appears that xo — $o- 
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FIG. 5: Logarithm of the stress s as a function of l/\ for 
HL results at temperatures both above (open symbols) and 
below (filled-in symbols) the glass transition. The curves are 
the predictions of the STZ theory at the same temperatures. 

At the outset, it might seem plausible to set 9$ = xo 
and thereby assume that there is only one energy scale 
that determines both the thermal and effective glass tem- 
peratures. This simplification cannot always be correct, 
however, because we know from the ubiquitous appear- 
ance of transient stress peaks in stress-strain curves mea- 
sured at high strain rates that slow thermal quenches can 
bring glasses into states where 8q < x < Xo- The con- 
ventional interpretation of these stress peaks is that the 
as-quenched state of the system has a low Active tem- 
perature x with a correspondingly small population of 
STZ's, and that plastic flow cannot begin until enough 
deformation has occurred to increase x to a value of order 
Xo- The Active temperature cannot be below 9o, so stress 
peaks could not be produced by this process if 9 = xo- 
In this situation, therefore, we expect that the HL model 
would not exhibit stress peaks except, possibly, at very 
high strain rates where xil) 1S substantially larger than 
Xo- 

In exploring fits to the viscosity data, however, we find 
one intriguingly different possibility. We can fit the HL 
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viscosity data with p = 2 in Eq. (|2.17[) - the prediction 
of the excitation-chain theory - but this fit looks best 
if we reduce 9$ to about 0.1, i.e. half the previous esti- 
mate. This lowered transition temperature would mean, 
for example, that the HL curves of stress versus strain 
rate shown in Fig [6] for values of 9 down to 0.1, would 
bend over and exhibit viscous behavior at strain rates 
much smaller than those measured by HL. If that were 
the case, then xo would be substantially higher than 9o, 
and the transient behavior of the HL model might be 
more interesting than predicted at the end of the last 
paragraph. 

B. Comparisons over the Full Range of the HL 
Data 

In order to extend our analysis to temperatures above 
the glass transition (Kauzmann) temperature 9q and to 
arbitrary strain rates, we must reexamine our assump- 
tion that the effective Kauzmann temperature xo is in- 
dependent of the bath temperature 9. To see what is 
happening here, we take advantage of our simplifying as- 
sumption that (3 — 1 in Eq. (|2.4p and solve this equation 
for x- 

1 _ T(s,9)/x(q) + K P (9)/9 

X T(s,9) + K p(9) • { - } 

Note that, for 9 > 9q, we can take the limits q — > and 
9 — > 9q in different orders and get different answers. If 
we take the limit of vanishing q first, so that r(s, 9) van- 
ishes, then x - > However, if we let 9 go to 9q first, so 
that p{9) vanishes, then x ~ * X.(q) ~* Xo f° r sufficiently 
small q. We see no a priori reason why the crossover 
between these limiting behaviors described by Eq. (|3.8j) 
cannot be physically correct for a temperature indepen- 
dent xo; but that assumption is qualitatively inconsistent 
with the HL data. Equation (|3.8|) implies that x always 
must lie between 9 and x(<z); but we have found (in a 
calculation not shown here) that, if we choose xo = $o to 
be a ^-independent constant, then the HL data in Fig[2] 
for 9 = 0.31 and 0.44 lie above the allowed range in the 
transition region between the small-g and large-g limits. 

A physically plausible alternative is to assume first that 
Xo > $o (the inequality being consistent with the exis- 
tence of transient stress peaks), and then that x(q ~^ 
0) = 9 whenever 9 exceeds xo- That is, 

lim *(<?) = / *° !^ <X0 ' (3.9) 
q ^o AyH ' [9 if 9 > xo- 

With this assumption, x(q ~ * 0) is a continuous function 
of 9 across the glass transition and, for 9 > xo, there 
are no small-g steady-state solutions other than x = 9. 
Nothing in this hypothesis precludes xo itself from being 
a function of temperature. In fact, we think we see a hint 
in the HL data that xo may be a decreasing function of 9 
up to the point where 9 = xo> beyond which x{q ~ * 0) = 
9. 



The assumption in Eq. (|3.9p accurately fits the HL 
data. In implementing Eq. (|3.9[) . we have simply rescaled 
the parameters xa and xi i n Eq. (|2.18[) so that the ra- 
tios xa/xo and Xi/Xo remain the same in both low and 
high temperature regions. We also have chosen n = 5 in 
Eq.flUJ and fl = 0.3 in Eq. (|2.12|) . Our results are shown 
in Figs. 1(3 and [3 

Figure [2] shows x as a function of strain rate q for a 
set of temperatures 9 both above and below 9q. As ex- 
pected, the HL data and STZ curves collapse to roughly 
a single curve for temperatures below 9q- Above 9q, as 
predicted by Eq. (|3.8[) , x —> & in the limit of small q; and 
the crossover to large- q behavior is correctly predicted by 
the theory. Figure [3] shows the same data as in Fig. [5] 
but here as the Arrhenius plot, Fig. [TJd, with the high- 
temperature data included. We show this different view 
of the data to emphasize that even the high-temperature 
values of x collapse to a single Arrhenius line, extrapo- 
lating to a single "melting" strain rate qo, in the limit of 
infinite x- 

Similarly, the graphs of stress s as functions of 1/x 
in Fig. [5l and those for s as functions of strain rate in 
Figs. Eta) and (b), show good agreement between theory 
and simulations at all temperatures and stresses. As in 
the preceding figures, these curves and the corresponding 
HL data points are shown at a selection of temperatures 
above and below 9q. Note, in Fig. [51 that the curves 
for 9 > 9q make smooth transitions to linear viscous be- 
havior at small strain rates. As expected, the curves for 
9 < 9q level off at temperature-dependent yield stresses. 
The short-dashed lines in Fig. EJa) illustrate the bifurca- 
tion at #o- At temperatures just above 9q, the material 
appears to exhibit a yield stress down to very small strain 
rates, but eventually flows even in the limit of arbitrar- 
ily small stress, thus showing how hard it is to make an 
accurate estimate of 9q by going to smaller and smaller 
strain rates. 

Two material-specific parameters, fl and s±, are best 
constrained using the HL data for \og(s/q) vs. log(g), 
shown in Fig. [4] The small strain-rate asymptote of 
these curves for temperatures above the glass transition 
is log(r)(9)), and the way in which the stress crosses over 
from its power-law dependence at large q to the yield 
stress or Newtonian viscosity at small q is strongly sen- 
sitive to STZ transition rate parameters fl and s±. We 
find that fl = 0.3 and Si = 0.08 fit the crossover behavior 
quite well. 

The curves in Fig. Q] show the STZ solutions for 
log(s/g) vs. log(g) at temperatures above and below the 
glass transition. For temperatures above the glass tran- 
sition, f](q) — s/q is very similar to the Cross form used 
by HL to extrapolate their data. However, the STZ solu- 
tions for r)(g) asymptote to a slightly smaller value of r\ 
as q — > than the one specified by HL. It is unclear how 
to best extrapolate this stress/strain rate data to s = 0, 
and this uncertainty explains why the HL viscosity values 
in Fig. 0Jl are slightly larger than the STZ fit. 
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FIG. 6: (a) Theoretical log-log plots of stress s as a function 
of strain rate q for a range of temperatures above and below 
the glass transition. The short-dashed curves are for 8 = 0.28 
and 6 = 0.30, both just above the glass transition temperature 
6*o = 0.2. (b) The same as (a) but for a smaller range of strain 
rates and including HL data at the selected temperatures 



IV. CONCLUDING REMARKS 

The natural way in which the STZ theory accounts 
for the wealth of data provided by the Haxton-Liu sim- 
ulations lends credibility to major elements of that the- 
ory, especially the form of the expression for the plastic 
strain rate in Eq. (|2.ip and the central role played by the 
effective disorder temperature. The self-consistent eval- 
uation of parameters in Section Mil depends directly on 
the first of these elements; and the remarkable appear- 
ance of an effective glass transition in the nonequilibrum 
glassy state seems to confirm the second. Nevertheless, 
unanswered questions abound. We conclude by listing 
some of them. 

How close is the analogy between the thermal and 
effective-thermal glass transitions? The analogy does 
seem surprisingly close, but there are interesting and po- 
tentially very important differences. We think that the 
near equality of the transition temperatures 8q and xo, if 
correct at all, must be an artifact of the HL model. We 
have argued here that 9q must be substantially smaller 
than xo m many realistic circumstances; and we suspect 
that this inequality might even turn out to be true in the 
HL model if measurements could be made at apprecia- 
bly smaller strain rates. It already is apparent that the 
effective Arrhenius activation energy A ks Tz in Eq. (|1.2p 



is larger than the thermal activation energy ks Tz by a 
factor of about 1.5; thus there is no universality in the 
energy scales. Why not? 

We see an even deeper question in this regard. As 
discussed in the paragraph preceding Eq. (|1.2|) , the fact 
that x{l) goes to a nonzero value xo in the limit of small 
q seems to have a robust explanation in purely kinetic 
terms. Is a similar argument relevant to the thermal 
glass transition? Would it imply that the viscosity - the 
analog of g _1 - rigorously does diverge at #o? 

What should the STZ transition rate R(s) really look 
like? Which features of it are model dependent? Which 
have some degree of universality? In Eq. lj2.12p . we have 
tried to go part of the way toward deducing R(s) from 
first-principles, many-body physics. Doing this correctly 
would be analogous to using statistical mechanics to com- 
pute the viscosity in the Navier Stokes equations. The 
problem here is harder, of course, because we do not yet 
have a precise definition of an STZ. Nevertheless, this 
problem needs to be given more thought in view of the 
new information we have obtained from the HL results. 

One especially intriguing aspect of this question 
emerges from our observation that R(s) exhibits Bag- 
nold scaling at large s. Under what circumstances may 
we expect this to happen? Is it special to the HL model? 
Ordinarily, we think of Bagnold scaling as being relevant 
to the large-stress behavior of the strain rate as a whole, 
i.e. to q(s) in present notation; whereas our q(s) be- 
haves quite differently. The special feature here is that 
the Bagnold behavior is modified by the STZ density, 
i.e. the factor exp(— 1/x) in Eq. (|2.ip . which is strongly 
q dependent. We suspect that the full Bagnold behav- 
ior would occur at strain rates larger than qo 7 where the 
STZ density is not a meaningful quantity, but where the 
HL model might look like a rapidly flowing granular ma- 
terial. A better understanding of this large-g situation 
might provide insight regarding the applicability of STZ 
theory to granular systems. 

Will the present version of the STZ theory accurately 
predict the time- dependent transient behavior of the HL 
model? One of the principal successes of earlier appli- 
cations of the STZ theory has been that it correctly ac- 
counts for response times observed, for example, in tran- 
sient stress-strain experiments. The key idea here is that 
there is a separation of time scales between slow pro- 
cesses, i.e. relaxation of plastic flow or disorder temper- 
ature, as opposed to fast processes like the response to 
perturbations of the STZ orientation variable m that ap- 
pears in Eq. (|2.8p . That separation seems problematic at 
the upper end of the scale of strain rates used by HL. We 
have no problem with it here because we consider only 
steady-state behavior. 

It will be interesting to look at transient behavior in 
the HL model to see where the fully dynamic version 
of the STZ theory may break down. Before doing that, 
however, it will be important to look at less esoteric be- 
havior, e.g. transient responses to changes in driving 
forces of systems initially quenched at different rates to 
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temperatures near the glass transition. Do such systems 
exhibit stress peaks? The model that we have used for 
interpreting the HL data, in which xo = implies that 
they should not occur here. Might transient experiments 
reveal that, as we suspect, xo > ^o? 

What are the implications for strain localization? The 
effective temperature provides a mechanism for strain lo- 
calization in amorphous materials - regions with a higher 
effective temperature are more likely to undergo plastic 
deformation, which generates heat that increases the lo- 
cal effective temperature. Although previous studies as- 
sumed a constant value for x UM, Eq. (|f .2[) shows that 
this assumption is only valid near the effective glass tran- 
sition. The fact that x increases with increasing strain 
rate is likely to enhance localization and may permit 
steady state localized solutions. 

What role do thermal fluctuations play in determin- 
ing the properties of the effective glass transition? We 
already have concluded that x{l ~ y 0) must be equal 



to 6 at high temperatures, 9 > xo', but we have lit- 
tle knowledge about what happens to the corresponding 
Xo when 8 < xo, and have simply assumed that it is a 
constant. More detailed information about this quantity 
might help us understand the molecular mechanisms that 
are at work here. 
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